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Abstract 

We prove that the number of monomer-dimer tilings of an n x n square 
grid, with m < n monomers in which no four tiles meet at any point is 
m2 m + (m + l)2 m+1 , when m and n have the same parity. In addition, we 
present a new proof of the result that there are n2 n ~ l such tilings with n 
monomers, which divides the tilings into n classes of size 2 n_1 . The sum of 
these tilings over all monomer counts has the closed form 2 n-1 (3n — 4) + 2 
and, curiously, this is equal to the sum of the squares of all parts in all 
compositions of n. We also describe two algorithms and a Gray code ordering 
for generating the n2 n ~ 1 tilings with n monomers, which are both based on 
our new proof. 

Keywords: tatami, monomer-dimer tiling, combinatorial generation, Gray 
code, enumeration, polyomino 
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1. Introduction 

Tatami mats are a traditional Japanese floor covering, whose outside is 
made of soft woven rush straw, and whose core is stuffed with rice straw. 
They have been in use by Japanese aristocracy since the 12th century, and 
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are so integral to Japanese culture, that they are often used to describe the 
square footage of a room. A standard full mat measures 6' x 3', and a half 
mat is 3' x 3'. 

An arrangement of the mats in which no four meet at any point is often 
preferred. We call such arrangements monomer-dimer tatami tilings. The 
monomer-dimer tiling in Figure [T] violates the tatami condition, and the one 
in Figure [2] does not. 

Figure 1: The tatami condition is that no four tiles may meet at a point. All such violations 
are circled. 



Figure 2: A tatami tiling of the 10 x 31 grid with 10 monomers, 64 vertical dimcrs and 
86 horizontal dimcrs. Compare this with Figure [5] 

Tilings with polyominoes are well studied, and they appear in physical 
models, the theory of regular languages, and combinatorics (see 0, [l?], E, H 



It is natural to ask, given a room's size, how many different arrangements 
of mats are there? Motivated by an exercise posed by Knuth ([lo|), Ruskey 
and Woodcock published the ordinary generating functions for fixed height, 



dimer-only tatami tilings in [14 



In Fall 2009 our research group discovered the structure which is published 
in and described in Section [1.11 Problems such as the enumeration of 
tatami tilings with r rows, c columns, and m monomers, have since become 
within reach. This paper completes such an enumeration for square grids. 

Alhazov et al. followed up on the aforementioned dimer-only research 
with a treatment of odd-area tatami tilings which include a single monomer. 
They closed [l[ with this remark: 
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However, the variety of tilings with arbitrary number of monomi- 
noes is quite "wild" in sense that such tilings cannot be easily 
decomposed, see Figure 11; therefore, most results presented here 
do not generalize to arbitrary number of monominoes, the tech- 
niques used here are not applicable, and it is expected that any 
characterization or enumeration of them would be much more 
complicated. 

The structure we found, however, reveals the opposite; the tilings with an 
arbitrary number of monomers are easily decomposed. The decomposition 
has a satisfying symmetry, it is ammenable to inductive arguments, and it 
shows that the complexity of a tatami tiling is linear in the dimensions of the 
grid (compare Figure [2] with Figure E]). We use it extensively to prove our 
main result. Let T(n, m) be the number ofnxn tatami tilings with exactly 
m monomers. 

Theorem 1. If n and m have the same parity, and m < n, then T(n,m) = 
m2 m + (m + l)2 m+1 . 

This confirms Conjecture 3 in (i[ for d = (as was promised there), and 
completes the enumeration of tatami tilings of square grids, since T(n, m) = 
when m > n. 

The proof that T(n,n) = n2 n_1 , in jij, counts tilings recursively. Due 
to the symmetry of the solution it is natural to search for a way to divide 
the tilings into n classes of size 2 n_1 ; indeed, Knuth mentioned this in a 
correspondence with Ruskey (jsj). We remedy the lack of symmetry in the 
original proof by presenting a new one which includes the desired partition 
and counts the tilings directly, rather than recursively. 

Nice round numbers such as these tend to be easily identified in out- 
put generated by exhaustively listing all tilings on a given size of grid. In 
a personal communication Ruskey (pj) Knuth noted there appeared to be 
2 n ~ 1 (3n — 4) + 2 tatami tilings of the n x n grid. We prove this result, and 
show that it is equal to the sum of the squares of all parts in all compositions 
of n. 

In Section [3] we use the partition in Theorem [2] to exhaustively generate 
all n x n tilings with n monomers, for even n, in a constant amount of time 
per tiling. The problem of generating the tilings in each of these partitions 
reduces to generating all bit strings of a fixed length, so our first algorithm 
is a straightforward application of this. The second algorithm generates the 
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tilings in a Gray code order, where the Gray code operation is the diagonal 
flip (defined in Section ITTTj) . This is acheived by generating the partitions in 
a certain order, and using a binary reflected Gray code to generate the bit 
strings. 

1.1. Structure 

Tatami tilings of rectangular grids have an underlying structure, called 
the T-diagram, which is an arrangement of four possible types of features, 
shown in Figures [3JHJ up to rotation and reflection (see |5| for a complete 
explanation). Each feature consists of a source, shown in colour, which forces 
the placement of up to four rays. Rays propagate to the boundary of the 
grid and do not intersect each other. 

A valid placement in the grid of a non-zero number of features determines 
a tiling uniquely. Otherwise, the trivial T-diagram corresponds to the four 
possible running bond patterns (brick laying pattern). 

(a) A loner (b) A vee feature, 

feature. 

Figure 3: These two types of features must have their coloured tiles on a boundary, as 
shown, up to rotation and reflection. 




(a) A vertical (b) A counter- 

bidimer feature. clockwise vortex 

feature. 



Figure 4: These two types of features may appear anywhere in a tiling provided that the 
coloured tiles are within the boundaries of the grid. 
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Monomers occur on the boundary of the grid, and inside vortices, but 
nowhere else. The T-diagram is also a partition of the tiles into blocks of 
horizontal and vertical running bond. The tiling in Figure |5] uses all four 
types of features, displays this partition into horizontal and vertical running 
bond. Its T-diagram is shown on its own in Figure [6j 

r~i i n i 1 

] 



i_i i_i i i i 

Figure 5: A tiling showing all four types of sources. Coloured in magenta, from left to 
right they are, a clockwise vortex, a vertical bidimer, a loner, a vee, and another loner. 




Figure 6: A T-diagram 



Throughout this paper we consider tilings on the nx n integer grid, with 
the origin at the bottom left of the tiling. Let the coordinate of a grid square 
be the point at its bottom left corner as well. 

A diagonal is a rotation of the following: a monomer at (xq,0), and 
(vertical) dimers covering the pairs of grid squares (xo + 1 + k,k) and (x + 
1 + k, k + 1), for each non-negative k such that xo + l + k < n — 1. A diagonal 
can be flipped in place, so that the monomer moves to (n — 1 , n — xq — 1 ) , and 
the dimers change orientation. This operation preserves the tatami condition. 
Diagonals are used and illustrated frequently in the next section. 

2. Square grids 

Let T n m be the set of n x n tatami tilings with m monomers. To find 
|T n)m |, we use a setup similar to that of Theorem 2 in jEJ. These lemmas 
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on the composition of such tilings seem apparent from Figures EHH1 but we 
prove them in general. 

We show that any tiling in T n>m with m < n has exactly one bidimer or 
vortex and we give a one-to-one correspondence between m and the shortest 
distance from this source to the boundary. Such a feature determines all tiles 
in the tiling except a number of diagonals that can be flipped independently. 
Proving the result becomes a matter of counting the number of allowable 
positions for the bidimer or vortex, each of which contributes a power of 2 
to the sum. 

For example, the 20 x 20 tiling in Figure [7] has a vertical bidimer which 
forces the placement of the green and blue tiles, while the remaining diagonals 
are coloured in alternating grey and magenta. There are eight such diagonals, 
so there are 2 s tilings of the 20 x 20 grid, with a vertical bidimer in the position 
shown. Each of these 2 8 tilings has exactly 10 monomers. 




bl x f -y f y f y f n-Xf- y f br 



Figure 7: Vertical bidimer. 

By Corollary 1 in j^j, a tiling in T n , jm , with m < n, must have at least 
one bidimer or vortex, which we will call /. Let (xf,yf) be the center of / in 
the cartesian plane, where the bottom left of the nx n tiling is at the origin. 

Define Xf as the lines (x — Xf) + {y — yf) = and (x — xf) — {y — yf) = 0, 
bounded by the grid boundary, which form an X through (xf, yf) (see Figures 
GH]). 
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bl x f -y s y f y f n-x f - y f br 



Figure 8: Horizontal bidimcr. 

Without loss of generality, we may re-orient a tiling, by rotating and 
reflecting, so that yf < x f < n/2. The upper arms of Xf intersect the left 
and right boundaries of the grid, while the lower arms intersect the bottom. 
In this range, the distance from this feature to the boundary of the grid is 

Vf- 

Lemma 1. Let T e T n ,m, with m < n. Then T has exactly one bidimer or 
vortex, f , but not both, and no vees. 

Proof. It is easy to see that if a is another bidimer or vortex in the tiling, 
then X a intersects Xf, which contradicts the fact that rays do not cross. A 
vee has the same rays as a bidimer, by replacing the adjacent monomers with 
a dimer, so the only features that may appear, besides /, are loners. □ 
Let Tf be the tiling with f as its only feature. This corresponds to the 
featurless "trivial" tilings in [5j. 

Lemma 2. The tiling Tf can be obtained from any tiling containing the 
feature f via a finite sequence of diagonal flips in which each monomer moves 
at most once. Reversing this sequence gives the original tiling. 

Proof. This is a simple modification of the proof of Lemma 3 in \^ . □ 
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bl x f -y s y f y s n-x f - y f br 



Figure 9: Counter clockwise vortex. Note that xj and yj are not integers. 

Lemma 3. All the diagonals in Tf can be flipped independently. 

Proof. We show that no two diagonals intersect, by sharing a monomer, or 
other grid squares. 

Let a and (3 be the monomers in distinct diagonals, possibly with a = (3, 
which intersect. Let L a be the line through a, and let Lp be the line through 
0. Their slopes must differ, because in sharing grid squares, they must share 
whole tiles. Thus we may assume that L a has slope 1, and Lp has slope — 1. 
Then L a and Lp intersect inside the grid, and therefore at least one of them 
instersects with Xf, which makes the corresponding diagonal unflippable. 
Therefore distinct diagonals may not intersect. □ 

The crux of the argument in Theorem [1] is this: 

Lemma 4. If Tf G T n ,m, with yf < Xf < n/2, then m = n — 2yj if f is 
a bidimer, and m = n — yf + 1 if f is a vortex. The number of (flippable) 
diagonals in Tf is 

n-2y f -2, ify f <x f ; 
n — 2yj — 1, if yf — Xf < n/2; and 
0, otherwise. 
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Proof. Let tl, tr, bl, and br be the segments of the boundary of the grid 
delimited by Xf, as shown in grey in Figure 10(a) We count the number of 



monomers and diagonals in Tf by measuring the lengths of tl, tr, bl, and br. 



minus the tiles in the rays of /, casified in Figures 10(b) -10(e 



n-Xf-fy 


Itl 


tr I 


Xf-Vf 




\ / 










br 





x f -y f n-Xf-fy 
(a) The range yf<Xf<n/2is green. 
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(b) 
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(c) 



(d) 



(e) 



Figure 10: The segments tl, tr, bl, and br. When tl and br have positive length, one of 
their combined tiles is part of a ray of /, shown in green or blue. Similarly for tr and bl 
combined. 



Figures 10(b) 10(e) show (in general) that if tl and br are non-zero, then 
together they contain exactly one grid square covered by the ray of /. The 
same is true of tr and bl. 

Because of the above, non-zero tl and br differ in parity, so their sum is 
odd, and therefore they must contain exactly one corner monomer (which is 
not in a diagonal). Similarly for tr and bl. 

We tabulate the numbers of monomers and diagonals in tl, tr, bl and br 
for different (xf,yf) in Tabled! and add pairs of rows to prove each case of 
the theorem statement. 

If yf < Xf, then m = Xf — yf + n — Xf — yj = n — 2yf if / is a bidimer, 



and m = n — 2?/ j + 1 if / is a vortex. The number of diagonals is n — 2yj — 2, 
as required for this case 

If yf = Xf < n/2, then m = + n — yf — yf = n — 2yj if / is a bidimer, 
and n — 2yf + 1 if / is a vortex. The number of diagonals is n — 2yf — 1. 
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(xf,yf) Positions Monomers Diagonals 



Uf = Xf bl and tr 

Uf < Xf bl and tr 

Uf = Xf = n/2 tl and br 

Uf + Xf < n tl and br 



n — Xf — yj n — Xf — yf — 1 




Table 1: 



If Uf = Xf = n/2, then m = if / is a bidimer, and m = 1 if it is a 



We sum up the number of positions which give a particular m in Lemma 
IU to count all such tilings. 

Theorem 1. If n and m have the same parity, and m < n, then T(n,m) = 
m2 m + (m + l)2 m+1 . 

Proof. We count the number of tilings Tf in T njm , for each bidimer or vortex, 
/, satisfying the conditions of Lemma HJ and for each of these we count the 
tilings obtainable via a set of (independent) diagonal flips. 

If / is a bidimer, then m = n — 2k, where k is the shortest distance from 
(xf,yf) to the boundary of the grid. These are the positions 

(k, k), (k, k + 1), {k, k + 2), . . . (A;, n - k), (k + 1, n - k), . . . , 

{n — k, n — k), {n — k, n — k — 1), . . . {n — k, k), {n — k — 1, k), . . . {k + 1, k), 

of which there are 4(n — k — k). We apply Lemma S] by making the necessary 
rotations and reflections, so that when m > 0, the four positions (k, k), (k, n— 
k),(n — k,n — k) and (n — k,k), have m — 1 diagonals, while the remaining 
4m — 4 positions have m — 2 diagonals. 

The same logic applies when / is a vortex and m > 1, and the results are 
in Table H 

By Lemmas [2] and [31 all of the tilings which contain / can be obtained 
from Tf by making a set of independent flips. Thus, the number of these is 
2 d (/) ; where d(f) is the number of diagonals in Tf. 

Lemma [1] tells us that there is no other way to obtain a tiling in T n>m , 
so we conclude by summing the 2 d ^s for each / such that T f G T n>m . Each 
term in the following sum comes from a row of Tabled in the same respective 
order, and similarly, the three factors in each sum term come from the last 



vortex. The number of diagonals is 0. 



□ 
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Type 



Feature 



Positions Diagonals 



x j = y.f 

x f = Vf 
x f < Vf 
x f < Vf 



h and v bidimers 
cc and c vortices 
h and v bidimers 
cc and c vortices 



4 m — 1 

4 m — 2 

4(m — 1) m — 2 

4(m — 2) m — 3 



Table 2: Horizontal, vertical, counterclockwise and clockwise are abbreviated as h, v, cc 
and c, respectively. We assume m > if / is a bidimcr, and m > 1 if / is a vortex. 

three columns of the table. 

T(n, m) =2 ■ 4 ■ 2 m ~ x + 2 • 4 • 2 m ~ 2 + 2 • 4(m - 1) • 2 m ~ 2 + 2 ■ 4(m - 2) • 2 m " 3 
=2 . 2 m+1 + 2 • 2 m + (m - l)2 m+1 + (m - 2)2 m 
=m2 m + (m + l)2 m+1 . 

That is, there are m2 m tilings with vortices and (m + l)2' m+1 with bidimers. 
This also satisfies T(n, 0) when / is a bidimer, and T(n, 1) when / is a 
vortex. □ 

This completes the enumeration ofnxn tatami tilings with m monomers. 

We present an alternate proof of Theorem 2 in [^j . The new proof orga- 
nizes the n x n tatami tilings with n monomers into n classes of size 2™" 1 , 
answering a challenge of Knuth, as noted in {(J . In addition to showing a sat- 
isfying symmtery, the resulting partition is used in Section [3] for exhaustively 
generating all such tilings in constant amortized time. 

The proof in j^j uses induction on the dimensions of the grid to show that 
n2 n ~ 1 = 4S(n), where 



The recurrence arises from relating the conflicts in flipping pairs of diagonals 
in an (n — 2) x [n — 2) tiling with those of an n x n tiling. 

We reuse the setup in [5|. Every tatami tiling of the n x n grid with n 
monomers can be obtained by performing a sequence of flips (of diagonals) 
on a running bond, in which no monomer is moved more than once, and the 
corner monomers remain fixed. As an immediate consequence, every tiling 
has exactly two corner monomers which are in adjacent corners, and as such, 
it has four distinct rotational symmetries. 

Let T n be the n x n tilings with n monomers, whose upper corners have 
monomers. It is sufficient to divide the n2 n ~ 3 tilings of T n , into n classes of 



S(n) = 2 n ~ 2 + AS(n - 2), where S(l) 



1 



and S(2) = 1. 
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size 2 n_3 , and then use the rotational symmetries of these to obtain the n 
classes of size 2 n ~ 1 of all the tilings, T n n . 

In this context, a flipped diagonal or monomer refers to one which was 
flipped from its initial position in the running bond, in the sequence of flips 
referred to above. Each monomer is in exactly two diagonals in the trivial 
tiling, and the monomer is flipped at most once, so its status can be described 
by a ternary symbol (see Figure [T5j) . 

Putting these ternary symbols in a list provides a representation for ele- 
ments of T n , which contains the status of each monomer as either unflipped, 
or flipped in a given direction. This gives the sequence of flips required to 
determine the tiling (See Figures [TT1 - [T2^ . 

For even n, let L = (lo,l\, . . . , ln-2 x ) be a ternary string, with l t rep- 
resenting the iih monomer from the bottom, on the left side of the trivial 
tiling. Assign ^ = if the monomer is not flipped, ij = — 1 if it is flipped 
down, and ij = 1 if it is flipped up. The downward diagonal of a monomer 
is denoted and the upward diagonal is l\ , regardless of whether they can, 
or have been flipped. Similarly, we define R = (rn, n rn-2 1 ). Their 
concatenation, L ■ R determines the tiling by describing the diagonal flips, 
up the left side, then up the right side. 



h 

h 

h-- 

h 

h 

h 

h 

h 




Figure 11: This 18 x 18 tiling is represented by (0,0,0,0,0, -1,0,1) ■ (0,0,0, 1,0,0,0,0). 
The position of a monomer in the trivial tiling reveals the number of dimcrs in its diagonals. 
The flipped diagonals arc blue, and the diagonals r{ and r\ are green, but they are 
unflipped. As expected, d\%(r\) + di$(r\) = 3 + 14 = 18 — 1. 
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If n is odd the representation is similar, noting that there are I ^f^J 
non-fixed monomers along the top and f 21 ^] along the bottom. We label 
from left to right, by To = (to, h, . . . , ii n=aj ) and Bo = (6 , £>i, • • • , &|~n^2~| ), 
respectively, and let = —1 if it is flipped on its left diagonal. Similarly, we 
define t~" as the rightward diagonal of tj, and so forth. 



^- <> <> 

\fr \\ VI to \A to to 

en 



A 3 



4> A A 



<9^r?" \/ <9^ 0^ l^O^O 



Figure 12: This 17 x 17 tiling is represented by (0,-1,0,0,0,0,0) • (0,0,-1,0,0,1,0,0). 
The flipped diagonals are green, and the two blue diagonals are coloured to show tX and 
t~X , which are the left and right diagonals of the same monomer. As expected, dn (t^) + 
dn{t^) = 10 + 6 = 17- 1. 

Let d n (a) be the number of dimers in the diagonal a, also called the length 
or size of the diagonal. It is a function of the index and direction of a: 



d n (a) = < 



2i + 


1. 


if 


a 


= if or a = rf or a = 


n — 2i — 


2. 


if 


a 


= l\ or a = r\ or a = 


2i + 


2. 


if 


a 


= t*~; and 


n — 2i — 


3. 


if 


a 





Conflicting pairs of diagonal are identified in a constant number of oper- 
ations by using the cases listed in Table [3j Conflicts occur only in cases 1(c) 
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and 2(b), and case 1(c) is always a conflict. In the case of 2(b), there is a 
conflict if and only if the sizes of the diagonals sum to at least n. 



A pair of monomers may be 

1. on the same side and flipped 

(a) in the same direction; 

(b) away from each other; or 

(c) toward each other; or 
on different sides and flipped 

(a) in different directions; or 

(b) in the same direction. 



2. 



Pair Conflict 2(b) iff 

7^ r { 

i ' 3 

7? r T 

I ' 3 



n<2(i + j] 
n>2{t + j] 
n<2(i + j) 
n>2(t + j] 



2 
4 
3 
5 



Table 3: Diagonal conflict cases. 

Let T n (a) C T n , where a is a diagonal (one of lj, lj, rj, rj, bj*, b^, tj*, or 
tf~), be defined as the collection of tilings in T n in which a is the longest 
flipped diagonal] for each flipped diagonal b, distinct from a, we have d n (b) < 
d n (a). 

Let T n (0) be the set of tilings in which no monomer is flipped on its 
longest diagonal. Note the distinction between a monomer flipped on its 
longest diagonal, and the longest flipped diagonal in the tiling (see Figure 



Initial position of monomer 



Figure 13: Each monomer in the tivial tiling is in two diagonals, a and a. Note that 
d n (a) + d n (a) = n — 1 

Here is the new proof, using the notation defined above. 
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Theorem 2(0). T(n,n) = n2 n ~ l 

(n classes of size 2 n ~ 1 ). Let D be the set of diagonals in the trivial tiling, 
and for each monomer in a diagonal a, we denote its other diagonal by a (see 
Figure [T3|) . Let A = {a G D : d n (a) > d n (a)} and let U denote a disjoint 
union. We prove that 



T n = T n (0)U ( M T n ( 



and, for a G A, 



. (2 n 3 ,2 n 2 ,n-2), if n is even; and 
T n a,T n 0,L4 ) = { ' ' ' .' (1) 

(2 ,2 ,7i — l), it n is odd, 



which gives the correct sizes, since |T n | = (n — 2)2 n 3 + 2 n = (n — l)2 n + 

2 n-3_ 

An available diagonal in T n (a) is a diagonal that may be flipped in a 
tiling of T n (a), subject to the constraints that a £ AU {0} and if a ^ 0, 
then d n (6) < d n (a), and 6 does not conflict with a. 

It suffices to show that the available diagonals are pairwise non-conflicting, 
and independent, and that for each T n (a), there are log 2 (|T n (a)|) of them. 
See Table @] for all the cases, and Figures [HHT5] for examples. 

For odd n, the diagonals a and a of the monomer on the middle column 
are equal in length. Thus is is not flipped at all in T n (0), giving the number 
|T n (0)| = 2 n ~ 3 , however, a, a G A, so \A\ = n - 1. 

For even n, no such symmetry exists, so we have |T n (0)| = 2 n_2 , and we 
must partition T n (0) into halves to obtain n classes of size 2 n_3 . □ 

Corollary 1. The number of n x n tatami tilings is 2 n ~ 1 (3n — 4) + 2. 



Proof. By Corollary 1 in [5|, we have that T(n, m) = when m > n. Let 

T ( n ) = Em>0 T ( ?2 ' m )' SO tllat 

\n/2\ 

T(n) = n2 n ~ 1 + ^ ((^ - 2i)2 n ~ 2i + (n - 2i + l)2 n ~ 2i+1 ) , 
i=i 

and notice that the sum simplifies to 

n-l 



T(n) = nT- 1 + ^iT 



=i 
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j must satisfy these constraints 


n odd 
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conflict 2(b) 
j < Sf 3 - i 


u n {Oj ) < a n \H ) 
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ronflirt 2(h) 
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symmetry with T n (t^) 




'* \ 2 / 




symmetry with T n (6^) 




T n (0) 


d n (^) < d n (tT*) 
^3< n r k 


d n (t]>) < d n {ty) 


d n {bf) < d n (b~>) 


dn(b?) < d n (bX) 

1 " \ j / J 


n even 
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T„(tf) 
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conflict 1(c) 
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conflict 2(b) 
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d n {r]) < d n (lj) 
=>- j > ^ — i 


T„(tf) 


conflict 1(c) 
=>j < i 


d n (l}) < d n (lj) 


dn(rj) < d n (lj) 


conflict 2(b) 






symmetry with T n (lj) 




Tn(rJ) 




symmetry with T n (/J) 




T n (0) 




dn(l]) < d n (l]) 


d n {rj) < d n {r]) 
^3<^ 


d n {r]) < d n (rj) 
^3>^ 



Table 4: Each T„(a) has exactly log 2 (|T„(a)|) independently flippable diagonals. See 
Figures [T4Hl~5l for examples. 



Now we use the fact that 2 k + 2 k+1 + ■■■ + 2 n ~ l = 2 n - 1 - 2 k + 1 = 2 n - 2 k 
to rearrange the sum. 

71-1 

T{n) =n2 n " 1 + J^i2 i 

i=l 

n-1 

=n2 n - 1 + (n-l)2"-^2 i 
=2 n ~ 1 (3n-4) + 2 

□ 

By cross referencing the result of Corollary [1] with the On-Line Encyclope- 
dia of Integer Sequences ((if), a surprising correspondence with integer com- 
positions of n was discovered (see sequence A027992). No proof is included 
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(a) T 17 (ttr). (b) Ti 7 (0). Note that tj" and 

are not available. 

Figure 14: Odd n. Available diagonals are marked by solid lines. 



in the sequence reference so we have included one here which illuminates a 
symmetry in this correspondence. 

We restate the result of Theorem [1] in terms of the distance of a bidimer 
or vortex from the boundary of the grid. 

Lemma 5. For 1 < k < n — 1, the number of tatami tilings of the nxn grid 
with a vortex or bidimer a distance (k+l)/2 from the boundary is (n — k)2 n ~ k . 

Proof. In the proof of Theorem [1] it was shown that the number of tilings of 
the nxn grid with m monomers and a vortex is m2 m , and with a bidimer it 
is (m + l)2 m+1 . Let d be the minimum distance from the center of a feature 
to the boundary of the grid. By Lemma H] in the case of a bidimer m = n — 2d 
and in the case of a vortex m = 2d + 1. Thus, in either case, the number of 
tilings written in terms of d is (n — 2d+l)2 n ~ 2d+1 or equivalently, (n — k)2 n ~ k 
where k — 2d — 1. □ 

Proposition 1. The number of n x n tatami tilings is equal to the sum of 
the squares of all parts in all compositions of n. 



Proof. Note that the number of red boxes in the diagram in |13[, see Figure 
[T6| is the sum of the squares of all parts in all compositions of 5. We consider 
in general this visual representation and show that for an integer composition 
of n the number of red boxes is equal to the number of tilings of the nxn 



17 




(a) T 18 (4). (b) T 18 (0). 



Figure 15: Even n. Available diagonals are marked by solid lines. 




Figure 16: A pictorial representation of the squares of all parts in all compositions of 5. 

grid. We proceed by partitioning these red boxes into n classes and showing 
that the cardinality of each class is equal to the number of tilings with a 
feature (or no features) a distance k/2 from the edge of the grid, for some 
k G {2,3,...,n}. 

For each of the n x n arrays in Figure [TH view the top left corner as 
having the (x, y) coordinates (0, 0) and the bottom right corner as having 
the coordinates (n, —n). Now consider a grouping of these red boxes into n 
classes, where a box is in class -X^for < k < n — 1, whenever its center 
lies on one of the lines y = —x ± k. The number of red boxes in class X Q is 
the number of red boxes along the line y = —x, which is n2 n ~ l since there 
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are 2 n ~ x compositions of n and n red boxes along each central diagonal. By 
Theorem [5] this is equal to the number of tilings of the n x n grid with the 
maximum number of monomers, or equivalently, no bidimers or vortices. 

For 1 < k < n — 1, if a summand of size % > k appears in a composition of 
n, then there are 2(i — k) red boxes along the lines y = ±k and none if i < k. 
In jij the authors show that the number of summands of size 1 < i < n 
over all compositions of n > 2 is given by (n + 3 — i)2^ n ~ 2 ~ i \ and trivially a 
summand of size n occurs once. Thus the number of red boxes in class Xk, 
for 1 < k < n — 1, is 

n-l 

2(n — k) + ( n + 3 ~ z )( z ~ fc)2 (n-1-i) . (2) 

i=fc+l 

Expression (|2J) simplifies to (n — k)2 n ~ k (details are omitted) which, by 
Lemma is the number of tilings of the n x n grid with a feature a distance 
from the boundary where 1 < k < n — 1. □ 

3. Combinatorial generation algorithms 

In this section we describe two combinatorial algorithms for generating 
all elements of T n , for even n. Each algorithm is a natural extension of 
our proof of Theorem [21 and we use the notation and ternary representation 
described there. 

The first combinatorial generation algorithm arises most naturally from 
our new proof of Theorem [5] (the proof in \^ leads to another, recursive 
algorithm). For each T n (a), Table H] can be used to determine the set of 
available diagonals. Several constant amortized time (CAT) algorithms exist 
for generating all subsets, so our algorithm is also CAT. 

3.1. Gray code 

The second algorithm uses a Gray code whose operation is the diagonal 
flip, which, in our representation, is equivalent to incrementing or decrement- 
ing a single ternary symbol. We use a binary reflected Gray code for each 
class of size 2 n_3 tilings, and order the classes so that the last tiling of one 
class is one flip away from first of the next class. If n = 2 (mod 4), this order 
begins with 

T n (/g); T n (rt_2_ i ), T n (/|), T n (rt_2_ 2 ), . . . , T n (/J), T n (r|) 
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where i is the largest integer such that d n {l\) > d n (lj), and j is the smallest 
integer such that d n (rj) > d n (A). The next class is T n (0), which is followed 
by 

Tn( r o)) T n (lt_ 2 ), T n (r|), T n (/t_2_ ), . . . T n (rJ_ 1 ), T n (^ +1 ). 

2 1 2 x 

For example, when n = 18 (See Figure [TTj) . this order is 

T 18 (4), T 18 (rJ), T 18 (4), T 18 (r£), T 18 (4), T 18 (r£), T 18 (4), T 18 (rJ), 

and then T 18 (0), continuing with 

T 18 (rJ), T 18 (4), T 18 (rl), T 18 (/J), T 18 (rJ), T 18 (4), T 18 (rJ), T 18 (4). 

The case where n = (mod 4) is similar. 

The binary reflected Gray code on n bits begins with the bitstring 
and ends with 10 n_1 . We assign to unflipped monomers and 1 to flipped 
monomers, since there is only one available flip for each flippable monomer. 
In this way, each T n (a) can be generated by one of the CAT algorithms in 

a. 

To generate two consecutive classes, T n (a) and T n (6), they must be sep- 
arated by exactly one flip. Notice that if b ^ and a / 0, then b is an 
available diagonal in T n (a), since a and b are on different sides, in different 
directions, and d n (b) = d n (a) — 1. In this case, we associate b with the 1 in 
the last bitstring, 10 n ~ 3 , of T n (a), and unflipping a allows us to begin gener- 
ating T n (6). If b = 0, then the order for the bitstring of T n (a) is arbitrary, 
however, T n (0) must be generated in reverse, ending with the trivial tiling. 
If a = 0, then T n (0) ends with the trivial tiling, so b is flipped and T n (b) 
is generated. 

The transitions from T n (a) to T n (6) are done in constant time, so this 
Gray code generation algorithm for T n runs in constant amortized time by 
making a sequence of calls to one of the CAT algorithms for generating binary 
reflected Gray codes in [jjj]. 

Two problems remain open in this discussion. Is there a Gray code for T n 
whose first and last elements also differ by 1 operation? A positive answer 
would provide a Gray code for T^ 11 . 

Our present discussion is about m = n and even n, but much of this 
applies to tilings of odd dimension or with less than n monomers. 
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4. Further research 



This paper answers some questions central to the enumeration and ex- 
haustive generation of square monomer- dimer tatami tilings, and similar find- 
ings are in preparation for the r x c grid, with r < c (see jij for related 
conjectures). The present research opens up the following paths: 

• find a useful representation for tilings in T n m which is ammenable to 
making diagonal flips; 

• exhaustively generate all of T n ,m in constant amortized time; 

• find the corresponding Gray code of Section 13.11 for odd n; 

• find an aesthetically pleasing bijection from the n x n tatami tilings 
to the red boxes of Figure [TBI f° r a U n - Preferably one which uses 
elements of the proof of Proposition [TJ and 

• explore n x n tilings with exactly v vertical dimers. 

This last point is in preparation, and might have been published here, 
save for its length. We will give a CAT algorithm for generating all of the 
n x n tatami tilings with exactly m monomers, v vertical dimers and h 
horizontal dimers, and we will explain the origin of the generating polynomial 
in Conjecture 4 in 
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